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\  ABSTRACT 

is  is  the  first  paper  of  a  series  in  which  we  outline  a  theory  of 
optimal  incomplete  block  designs  for  comparing  several  treatments  with  a 
control,  and  give  plans  of  optimal  designs  to  use  for  this  purpose.  In 
the  present  paper  we  develop  the  basic  theory  of  the  incomplete  block  de¬ 
signs  which  are  appropriate  for  the  many-to-one  comparisons  problem;  the 
optimality  considerations  for  such  designs  will  be  addressed  in  subsequent 
papers.  For  this  problem  we  propose  a  new  general  class  of  incomplete 
block  designs  which  are  balanced  with  respect  to  test  treatments.  We 
shall  use  the  abbreviation  BTIB  to  refer  to  such  designs.  We  study  their 
structure,  and  give  some  methods  of  construction.  A  procedure  for  making 
exact  joint  confidence  statements  for  our  multiple  comparisons  problem 
is  described.  Using  a  new  concept  of  admissibility  of  designs,  it  is  shown 
how  "inferior"  designs  can  be  eliminated  from  consideration .  Some  open 


problems  concerning  construction  of  BTIB  designs  are  posed. 


Key  words  and  phrases:  Multiple  comparisons  with  a  control,  incomplete 
block  designs,  balanced  treatment  incomplete  block  (BTIB)  designs,  admis¬ 
sible  designs,  optimal  designs,  applications  of  equicorrelated  multivariate 
normal  and  multivariate  Student's  t  distributions. 


1.  INTRODUCTION  AND  SUMMARY 

In  many  industrial,  agricultural  and  biological  experiments  it  is 
often  desired  to  compare  simultaneously  several  test  treatments  with  a 
control  treatment.  The  earliest  correct  work  on  this  problem  was  car¬ 
ried  out  by  Dunnett  [8],  [9],  Dunnett  [8]  also  posed  (but  did  not 
solve)  the  problem  of  optimally  allocating  experimental  units  to  control 
and  test  treatments  so  as  to  maximize  the  probability  associated  with 
the  joint  confidence  statement  concerning  the  many-to-one  comparisons 
between  the  mean  of  the  control  treatment  and  the  means  of  the  test  treat¬ 
ments.  Bechhofer  and  his  coworkers  solved  this  optimal  allocation  problem 
in  a  series  of  papers  [1],  [2],  [3]. 

In  all  of  the  aforementioned  papers  it  was  tacitly  assumed  that  a 
completely  randomized  (OR)  design  was  used.  However,  many  practical  sit¬ 
uations  may  require  the  blocking  of  treatments  in  order  to  cut  down  on 
bias  and  improve  the  precision  of  the  experiment.  If  the  block  size  is 
large  enough  to  accomodate  one  replication  of  all  of  the  test  treatments 
and  additional  control  treatments  as  well,  then  the  design  and  analysis  of 
replications  of  the  experiment  can  be  carried  out  using  the  optimal  alloca¬ 
tions  described  in  [1]  and  [2]  with  only  the  usual  modifications. 

We  shall  study  the  multiple  comparisons  problem  in  the  situation  which 
commonly  occurs  in  practice,  i.e.,  when  all  of  the  blocks  have  a  common 
size  but  the  block  size  is  less  than  the  total  number  of  treatments.  Robson 
[21]  pointed  out  that  Dunnett' s  procedure  can  be  extended  to  the  case  in 
which  a  balanced  incomplete  block  (BIB)  design  between  all  of  the  treatments 
(including  the  control  treatment)  is  used.  Cox  [7,  p.  238]  noted  that  BIB 
designs  are  perhaps  not  appropriate  for  the  multiple  comparisons  with  the 
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control  (MCC)  problem  because  of  the  special  role  played  by  the  control 
treatment.  He  suggested  a  design  which  employs  the  control  treatment 
an  equal  number  of  times  (once,  twice,  etc.)  in  each  block,  the  test 
treatments  forming  a  BIB  design  in  the  remaining  plots  of  the  blocks;  no 
analytical  details  were  given  for  this  t ropcsed  design.  Pesek  [20]  has 
given  analytical  details  for  a  special  case  of  Cox's  design  (i.e.,  the  con¬ 
trol  treatment  is  employed  once  in  each  block);  he  shows  that  this  design 
is  more  efficient  than  a  BIB  design  for  comparisons  with  a  control  but 
it  is  less  efficient  for  pairwise  comparisons  between  the  test  treatments. 
It  should  be  noted  that  even  Cox's  more  general  design  is  quite  restrict 
tive. 

In  our  series  of  papers  we  shall  propose  a  theory  of  optimal  incom¬ 
plete  block  designs  for  the  MCC  problem,  and  give  plans  of  optimal  designs 
which  would  be  useful  in  practice.  The  present  paper  deals  with  the  basic 
theory  underlying  our  designs;  the  problem  of  choosing  an  optimal  design 
from  a  set  of  contending  admissible  designs  is  addressed  in  the  subsequent 
papers  of  this  series. 

In  Section  2  we  give  the  linear  model  considered,  and  the  underlying 
assumptions  and  notation.  We  define  the  class  of  BTIB  designs  in  Section 
3,  and  find  the  necessary  and  sufficient  conditions  that  a  design  must 
satisfy  in  order  to  be  a  member  of  this  class.  Several  methods  of  construe 
tion  of  such  designs  are  given.  Expressions  for  the  best  linear  unbiased 
estimators  (BLUE's)  of  the  treatment  differences  of  interest  for  the  MCC 
problem  are  given  in  Section  4  along  with  the  analysis  of  variance  table 
associated  with  our  design;  it  is  shown  how  one  can  obtain  joint  confidence 
interval  estimates  of  the  differences  of  interest  using  existing  tables  of 
the  integral  of  the  multivariate  t  or  the  multivariate  normal  distribution 
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The  concept  of  the  admissibility  of  a  3TIB  design  is  introduced  in  Section 
5;  it  is  shown  there  how  this  criterion  can  be  used  to  rule  out  inferior 
designs.  In  Section  6  some  remarks  are  offered  concerning  the  content  of 
future  papers,  and  several  open  problems  are  mentioned.  Mathematical  de¬ 
tails  of  the  results  in  Section  3  and  4  are  given  in  the  Appendix. 


2.  PRELIMINARIES 

Let  the  treatments  be  indexed  by  0,1,..., p  with  0  denoting  the 
control  treatment  and  l,2,...,p  denoting  the  p  £  2  test  treatments. 

Let  k  <  p+1  denote  the  common  size  of  each  block,  and  let  b  denote  the 
number  of  blocks  available  for  experimentation.  Thus  N  =  kb  is  the  to¬ 
tal  number  of  experimental  units.  If  treatment  i  is  assigned  to  the  hth 
plot  of  the  jth  block  (O^i^p,  l<.h<J<,  lfjlb),  let  denote  the  cor¬ 

responding  random  variable;  we  assume  the  usual  additive  linear  model  (no 
treatment  *  block  interaction) 


Y.  =  u  +  a.  +  6 .  +  e.  .. 
iih  i  j  lih 


(2.1) 


with  V  a.  =  £  8.  =  0;  the  e.  are  assumed  to  be  i.i.d.  N(0,o‘) 

i=0  j=l  :  i;jh 

random  variables.  It  is  desired  to  make  an  exact  joint  confidence  state¬ 
ment  (employing  one-sided  or  two-sided  intervals)  concerning  the  p 
differences  a  -  a.  based  on  their  BLUE's  a.  -  a.  (l<i<n). 

0  l  0  l  — 


3.  CHOICE  OF  THE  CLASS  OF  DESIGNS 

3 . 1  BTI3  designs 

Since  it  is  desired  to  make  a  confidence  statement  which  applies  si¬ 
multaneously  to  all  of  the  p  differences  aQ  -  a^  (l<j.<p),  we  shall 


Ml*'1 


Jtek 
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regard  our  problem  as  being  symmetric  in  these  differences.  To  this  end 

A  ^  2  2 

we  consider  a  class  of  designs  for  which  Var{aQ-a.}  =  n  a  /N  ( l^i  <p ) 

A  A  A  A 

and  Corr{an-a.  ,  a  -a.  }  =  o  (i  ^i  ;  l^i. ,i-^p);  the  parameters  n 

U  —  -  U  X,a  x  x  xx 

X  X 

and  p  depend  on  the  design  employed.  We  shall  refer  to  such  designs 
as  3TI3  designs  since  they  are  balanced  with  respect  to  the  test  treatments. 
The  following  theorem  states  the  necessary  and  sufficient  conditions  that 
a  design  must  satisfy  in  order  that  it  be  a  3TIB  design.  The  proof  of 
this  theorem  is  given  in  the  Appendix;  in  the  process  we  also  derive 
expressions  for  VaHa^-ci^}  and  p.  These  quantities  play  a  crucial 
role  in  our  later  considerations. 


Theorem  3.1:  For  given  (p,k,b)  consider  a  design  with  the  incidence  ma¬ 
trix  {r..}  where  r. .  is  the  number  of  replications  of  the  ith  treatment 
^  V  b  — 

in  the  jth  block.  Let  A.  .  =  V  r.  .r.  .  denote  the  total  number  of 

1112  j=l  11^ 

times  that  the  i^th  treatment  appears  with  the  i„th  treatment  in  the  same 
block  over  the  whole  design  (i  /i_;  O^i  ,i  <d).  Then  the  necessary  and 

XX  X  X 

sufficient  conditions  for  a  design  to  be  3TI3  are: 


xoP  =  xo  (say) 


A  .  =  A,  (say) 

P-1.P  1 


(3.1) 


In  other  words,  each  test  treatment  must  appear  with  (i.e.,  in  the  same 
block  as)  the  control  treatment  the  same  total  number  of  times  (AQ)  over 
the  design,  and  each  test  treatment  must  appear  with  each  other  test  treat¬ 
ment  the  same  total  number  of  times  (A^)  over  the  design. 

An  example  of  a  3TIB  design  for  which  (p,k,b)  =  (4,3,10)  and 
Aq  =  4,  A  =  2  is  given  by 


Additional  examples  of  some  selected  BTIB  designs  are  given  in  (3.3) 


Expressions  for  Var{a  -a.}  and  Corr{u  -a 


Remark  3.1:  We  note  that  Theorem  3.1  places  no  restriction  on  r.  =  ) 

1  j=l 

(l^i^?) ,  the  number  of  replications  of  the  ith  test  treatment,  and  hence 


a  design  can  be  BTIB  without  the  r.  (l<i.<p)  being  equal.  Such  a  design 


for  which 


2  with  r 


At  this  point  we  indicate  several  methods  of  constructing  BTIB  designs 


As  a  starting  point  we  introduce  the  concept  of  a  minimal  BTIB  generator 


design  (more  simply  referred  to  as  a  generator  design).  For  given  (p,k) 


a  generator  design  is  a  BTIB  design  no  proper  subset  of  whose  blocks  forms  a 


BTIB  design;  here  at  least  one  of  X  ,A  >  0.  Thus 
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are  BTIB  designs  with  ( )  =  (1,0),  (0,1),  (1,1),  respectively; 
however,  only  DQ  and  are  generator  designs  (as  are  the  designs 

given  by  (3.2)  and  (3.3)).  Design  D  of  (3.4)  suggests  the  role  of 
generator  designs,  for  given  (p,k)  there  are  several  generator  de¬ 
signs;  e.g.,  for  each  p  ^  2,  k  =  2  there  are  exactly  two  generator 
designs.  (See  (3.5),  below.)  By  taking  unions  of  replications  of  these 
generator  designs,  at  least  one  of  which  has  A^  >  0,  we  obtain  an 
implementable  BTIB  design.  The  problem  of  determining  how  many  genera¬ 
tor  designs  exist  for  arbitrary  given  (p,k)  is  an  open  problem  (it  can 
be  shown  that  there  are  a  finite  number)  which  appears  to  be  very  formid¬ 
able;  we  comment  on  this  problem  in  Section  6,  and  consider  it  in 
Bechhofer  and  Tamnane  [5].  In  the  sequel  we  consider  only  implementable 
BTIB  designs;  we  also  assume  that  no  block  contains  only  one  of  the  p+1 
treatments  since  clearly  such  a  block  contributes  no  information  to  the 
estimation  of  the  a„  -  a.. 

0  l 

Suppose  that  for  given  (p,k)  there  are  n  generator  designs 

(l^iiji).  Let  A^  ^ be  the  design  parameters  associated  with 

D.,  and  let  b.  be  the  number  of  blocks  required  by  D.  (l^isn).  Then  a 
n 

BTIB  design  D  =  u  f.D.  obtained  by  forming  unions  of  f.  >_  0  replica- 

i=l  n  ,  .  .  n  ... 

r  (l)  r  vl) 

tions  of  has  the  design  parameters  XQ  =  \  f^AQ  ,  X^  =  L 

n  i=l  i=l 

and  requires  b  =  T  f.b.  blocks.  The  set  of  D.  with  f.  >  0  will 
. L,  i  l  l  l 

1=1 

be  referred  to  as  the  support  of  D. 

The  following  is  another  example  of  generator  designs  (generalizing 
(3.4)).  As  noted  above,  for  each  p  >  2,  k  =  2  there  are  exactly  two  gen¬ 
erator  designs.  They  are 


Co  0 

l1  2 


1 

3 


(3.5) 
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From  these  generator  designs,  implemer.table  BTIB  designs  of  the  type 

D  =  f,D_  u  f,D,  can  be  constructed  for  f„  >  1,  f,  >  0-  the  corre¬ 
ct  0  1  1  0  ~  1  — 

sconding  design  parameters  for  D  are  ,  \  =  fn  while 

b  =  f0p  +  f1?(p-l)/2. 

We  now  consider  several  methods  of  constructing  generator  designs 
for  k  ^  3;  this  list  is  not  exhaustive. 


Method  I_:  The  above  example  suggests  the  following  method  of  construct¬ 
ing  a  class  of  generator  designs:  For  given  (p,k),  a  generator  design 
D  will  have  m+1  plots  in  each  block  assigned  to  the  control  treat¬ 
ment:  the  p  test  treatments  are  assigned  to  the  remaining  k-m-1  plots 
of  the  b  blocks  (05mSk-2)  in  such  a  way  as  to  form  a  BIB  design. 

The  generator  design  D,  contains  no  control  treatments;  it  consists 

K  ~  x 

of  a  R3  (3IB)  design  between  the  p  test  treatments  if  p  =  k  (p>k). 

Thus  for  (?,k)  =  (3,3)  we  have  the  following  three  generator  designs 
in  this  class. 


Co  0  o'! 


foool 


fll 


Do  s<  1 


D,  =  <  0  0  0 

±  ' 


<  ! 


l2 


3  3 


(3.6) 


Method  II:  Starting  with  a  313  design  between  t  >  p  treatments  in 
b  blocks,  one  can  relabel  the  treatments  p+l,p+2, . . . ,t  by  zeros  to 
obtain  a  new  design  which  is  a  3TI3  design  with  possibly  an  additional 
block  or  blocks,  each  one  of  the  latter  containing  only  one  test  treat¬ 
ment  or  only  the  control  treatment.  After  deleting  all  of  these  one- 
treatment  blocks,  and  identifying  the  support  of  the  resulting  BTIB 


design,  one  obtains  the  desired  generator  design(s).  As  an  example  we 
consider  the  following  3IB  design  for  (t,k,b)  =  (7,3,7): 


r 

1  2  3  4  5 

\  7  1  2  3  4 

3  4  5  6  7 

w. 


(3.7) 


3y  replacing  the  sevens  by  zeros  one  obtains  a  generator  design  for 
(p,k,b)  =  (6,3,7)  with  A^  =  1,  X  =  1.  This  process  can  be  continued 
by  then  replacing  the  sixes  by  zeros  to  obtain  a  generator  design  for 
(p,k,b)  =  (5,3,7)  with  XQ  =  2,  =  1;  continuing,  one  can  then  replace 

the  fives  by  zeros  to  obtain  a  generator  design  for  (p,k,b)  =  (4,3,7) 
with  A  =  3,  A  =  1.  Finally,  replacing  the  fours  by  zeros  one  obtains 
the  union  of  the  two  generator  designs  DQ  and  D  of  (3.6)  with  a 
block  containing  all  zeros.  The  3TIB  designs  obtained  in  this  way  for 
k  =  3,  b  =  7  are 

?  =  6,  b  =  7  (XQ  =  1,  \1  =  1) 


r 

0  0  0  1 

(12  4  2 

3  6  5  4 


D  =  5,  b  =  7  ( X 


f 


0  0  0  0 


1  2 

5  3 

6  5 


3 

4 
6 


=  2,  X1  =  1) 

0  1  2>1 


(3.3a) 


1  1  3  4  0  2  3 


(3.3b) 
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p  =  4,  b  =  7  (X 


=  3,  Xx  =  1) 

-\ 

0  0  1 
0  0  2  \ 

2  4  4 


(3.3c) 


L  J 

Remark  3.2:  From  (3.8a)  it  is  clear  that  every  BIB  design  (or  Youden 
square  design)  involving  t  treatments  yields  a  BTIB  design  with 
p  =  t-1  test  treatments. 


Remark  3.3:  It  is  well-known  (see,  e.g. ,  John  [15])  that  unequal  repli¬ 
cate  designs  having  the  same  properties  as  BIB  designs  can  be  constructed. 
Such  designs  are  therefore  BTIB  designs.  An  example  of  such  a  generator 
design  for  (p,k,b)  =  (5,3,8)  with  XQ  =  0,  X^  =  2  is  given  by 


f 

111212341 
22335555 
34445555 


(3.9) 


Method  III:  Consider  a  group  divisible  partially  balanced  incomplete 
block  (GD-PBIB)  design  with  two  associate  classes  between  t  treatments 
in  blocks  of  size  k.  The  association  scheme  of  such  a  GD-PBIB  design 
can  be  represented  in  the  form  of  an  mxn  array  (with  mn  =  t).  Any 
two  treatments  in  the  same  row  of  the  array  are  first  associates  and 
those  in  different  rows  are  second  associates.  Suppose  that  m  >,  k;  one 
can  then  take  p  =  m  and  relabel  the  entries  in  n^  >  0  columns  of  the 
array  by  l,2,...,p  and  the  entries  in  the  remaining  n^  =  n-n^  >  0 


columns  by  zeros,  thus  obtaining  a  BTIB  design.  As  with  Method  II,  such  a 
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design  may  not  be  a  generator  design  and  may  contain  some  blocks  which 
must  be  deleted.  After  deleting  such  blocks  a  BTIB  design  is  obtained. 

By  identifying  the  support  of  this  resulting  design  the  desired  genera¬ 
tor  design(s)  are  obtained;  some  (or  all)  of  these  can  usually  be  obtained 
by  the  previous  two  methods.  If  n  >  k  one  can  take  p  =  n  and  then 
relabel  the  entries  in  m1  >  0  rows  of  the  array  by  l,2,...,p  and  the 
entries  in  m„  =  m-m  >  0  rows  by  zeros,  thus  obtaining  a  BTIB  design, 
possibly  with  blocks  which  must  be  deleted. 

To  illustrate  the  use  of  this  method  consider  the  GD-PBIB  design  # R20 
(for  k=3 ,  t=12,  m=4,  n=3 ,  b=36)  in  the  monograph  by  Bose,  Clatworthy  and 
Shrikhande  [6]  which  has  the  following  association  scheme: 

-\ 

5  9 

6  10 

> 

7  11 

8  12  J  . 

By  relabeling  the  treatments  5  through  8  by  1  through  4,  and  9  through 
12  by  zeros  one  obtains  the  union  of  a  BTIB  design  with  a  design  contain¬ 
ing  one  block  with  only  zeros.  After  deleting  that  block  the  support  of 
the  remaining  BTIB  design  consists  of  the  following: 

r 

0  0  0  0  0 

2  replications  of  /  1  1  1  2  2 

2  3  4  3  4 


o'! 
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r 

0  0  0 

2  replications  of  <  0  0  0 

12  3 

L 

r 

ill 

2  replications  of  (  2  2  3 

3  4  4 

r 

0  0  0 

1  replication  of  l  1  1  2 

2  4  4 


0 
0 
4 

~\ 

2 

3  r 

4 

J 

Oil 
0  2  3 
3  3  4 


(3.10b) 


(3.10c) 


( 3 . lOd ) 


Thus  the  designs  given  by  (3.10a)-(3.10d)  are  generator  designs  for 
p  =  4,  k  =  3;  the  first  three  designs  are  obtainable  by  Method  I  while 
the  fourth  is  not.  Using  the  definition  of  admissibility  of  a  design 
as  given  in  Section  5,  we  shall  see  that  this  fourth  design  with  b  =  7, 
=  2,  X^  =  2  is  admissible  relative  to  the  design  obtained  for 
(p,k,b)  =  (4,3,7)  by  Method  II. 


Method  IV:  Suppose  that  for  given  (p,k)  we  have  a  generator  design 
D,  with  >  0.  Then  a  new  generator  design  D„  for  the  same  (p,k) 
can  be  obtained  by  taking  a  "complement"  of  P,  in  the  following  way: 
Separate  the  blocks  of  in  different  sets  so  that  each  block  in  a 

given  set  has  zero  assigned  in  an  equal  number  of  plots  (0  times,  1  time, 
etc.).  For  example,  consider  the  design  (3.9c)  the  blocks  of  which  can 
be  separated  into  three  sets  as  follows: 
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3  3  4 


0  0  0 

0  0  0 
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For  each  set  of  D,  write  its  "complementary"  set  (with  zero  assigned 
in  the  same  number  of  plots)  so  that  the  union  of  that  set  with  its 
complementary  set  forms  a  generator  design:  if  r„  =  0  or  1  (l£iip) 
then  that  union  is  simply  a  generator  design  that  can  be  constructed 
by  Method  I.  These  complementary  sets  in  the  present  example  are 
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by  taking  their  union  we  obtain  the  generator  design  D?  given  by  (3.10d). 
The  b-values  for  and  its  complement  are  not  in  general  equal , 

although  in  the  present  example  they  are. 


4.  JOINT  CONFIDENCE  STATEMENTS 
4 . 1  Expressions  for  estimates 

We  first  give  the  expressions  for  the  BLUE  a„  -  a.  of  a  -  a. 

0  i  0  l 

(l^i^p)*  Let  T.  denote  the  sum  of  all  observations  obtained  with  the  ith 

treatment  (0<^i<p),  and  let  B.  denote  the  sum  of  all  observations  in  the  jth 

D 

block  (l4i£b).  Define  B*  =  T  r..B.  and  let  Q.  =  kT.  -  B*  (0<i<p). 

i  id  :  i  11-=^ 

Then 


1  *  W  V3!  ,,  .  % 


(4.1) 


where 


Source  of 
Variation 


Treatments 
(adjusted ) 


Sum  of  Squares 


Qo{(p-1)(VAi)+(p2+3)Vi}  i=i 0i 

k(p+l)2\ ‘(i0+pAi)  k(AQ+pX1) 


1  b  9  r2 

-  y  B2  -  — 
k  j£i  i  N 


(by  subtraction) 


The  expressions  in  the  table  are  derived  in  the  Appendix;  the  symbols 


G  and  H  are  also  defined  there.  We  note  that  if 


X  =  X  (say). 


then  SS^  ^  ,,  reduces  to  V  Q./kX(p+l),  i.e.,  the 

treatments(adnusted)  .  „  1 

1  =  0 

same  expression  as  for  a  BIB  design;  this  latter  expression  thus  holds 
for  any  completely  balanced  design  (such  as  (3.3)  which  is  not  a  BIB 
design) . 

4.2  Confidence  statements 

Joint  100 (1-a)  percent  confidence  intervals  for  the  -  a.  (l^iip) 
are  given  below. 

I.  One-sided  confidence  intervals 


a)  For  o  unknown: 


a-  -  a.  >  an  -  a.  -  t^°^  s  n/ /fT  (l<i<D) 
0  i^0  l  v,p,p  v  =  — 


(4.6) 


In  (4.6)  t^Clp  p  denotes  the  upper  equicoordinate  a  point  of  the  D- 

variate  eauicorrelated  central  t-distribution  with  common  correlation  p, 

and  with  d.f.  v  (as  defined  by  Dunnett  and  Sobel  [10]).  Krishnaiah  and 

Armitage  [17]  have  tabled  t^0^  to  three  significant  figures  for 

v ,  p  9  p 

p  =  1(1)10;  v  =  5(1)35;  a  =  0.05,0.01;  p  =  0.0(0. 1)0. 9. 

2 

b)  For  o  known: 


a„  -  a.  >  a„ 
0  i  =  0 


a . 

l 


on/'^N  (l<i<p) 

p,p  = 


(4.7) 


In  (4.7)  (=  t^a^  for  v  =  °°)  denotes  the  upper  equicoordinate 

p  ,p  v,p,p 

a  point  of  the  p-variate  equicorrelated  standard  normal  distribution 
with  common  correlation  p.  Gupta,  Nagel  and  Panchapakesan  [13]  have 


tabled  z^a  to  five  significant  figures  for  p  =  1(1)10(2)50: 
p,P 

a  =  0.250,  0.100,  0.050,  0.025,  0.010;  0  =  0.1(0. 1)0. 9,  0.125(0.125)0.875 
1/3,  2/3.  Two  other  relevant  references  are  Gupta  [12]  and  Miltcn  [19]. 

II :  Two-sided  confidence  intervals 
2 

a)  For  o  unknown: 


a  -  ot .  <•:  [a  -a.±t'^a^  s  n//Tj]  (l<_i<p) 
0  l  0  l  v ,p,p  v  ~=* 


(4.8) 


(a ) 

In  (4.8)  t'  satisfies 

v,p,p 


P{  1 1 .  i  <  t'(a)  (l<i<p)} 

l  —  v  ,p  ,p 


=  1  -  a 


where  (t  , ...,t  )  has  a  p-variate  equicorrelated  central  t  distribu- 
1  P 


tion  with  common  correlation  p,  and  with  d.f.  v  (Dunnett  and  Sobel 


(a ) 

[10]).  i.e.,  t’  is  the  upper  a  point  of  the  distribution  of 

v,p,p 


max(|t.|  (l^i^p)} .  Hahn  and  Hendrickson  [14]  have  tabled  t1  to  four 


significant  figures  for  p  =  1(1)6(2)12,15,20;  v  =  3(1)12,15(5)30,40,60; 
a  =  0.10,0.05,0.01,  and  p  =  0.0(0. 2)0. 4, 0.5.  See  also,  Krishnaiah  and 


Armitage  [18]  who  have  tabled  (t'^°^  )^  to  two  decimal  Places  for 

v  ,p,p 


p  =  1(1)10;  v  =  5(1)35;  a  =  0.05,0.01;  p  =  0.0(0. 1)0. 9. 


b )  For  a  known : 


a.  -  a.  e  [a„-a.  ±  z'^a^ari//N]  (l<i<p)  (4.9) 

0  i  0  i  p,p  — 

,(cO  (_  t,(o) 

«  \i 


In  (4.9)  z 


for  v 


®)  satisfies 
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P{ |Z. I  £  (l£i£p)}  =  1  -  a 

l  -  p,P  - 


where  has  a  p-variate  equicorrelated  standard  normal 


distribution  with  common  correlation  p,  i.e.,  z' 


is  the  upper 


a  point  of  the  distribution  of  max(|Z^|  (l<i<p)}.  Hahn  and 
Hendrickson's  tables  [14]  with  v  =  60  can  be  used  here  to  obtain  a 
conservative  approximation  to  v  = 


5.  THE  CLASS  OF  ADMISSIBLE  DESIGNS 
5 . 1  Optimal  and  admissible  designs 

We  have  shown  how  the  experimenter  can  make  joint  confidence 
statements  concerning  the  aQ  -  (l£i<p)  when  the  experiment  was 

conducted  using  a  BTIB  design.  However  we  have  not  proposed  a  ration¬ 
ale  for  choosing  one  design  from  a  set  of  contending  BTIB  designs. 
Toward  that  end  we  define  an  optimal  design  as  that  design  which  for 
given  (p,k,b)  maximizes  the  appropriate  confidence  coefficient  (for 
one-sided  or  two-sided  common  confidence  intervals)  subject  to  a  spec¬ 
ified  limit  on  the  "width"  of  the  intervals.  (One  might  also  define 
an  optimal  design  as  that  design  which  for  given  (p,k)  and  specified 
joint  confidence  coefficient  (for  one-sided  or  two-sided  confidence  in¬ 
tervals)  minimizes  b,  the  total  number  of  blocks  required  for  the 
experiment,  subject  to  a  specified  limit  on  the  common  "width"  of  the 

intervals . )  The  solution  to  both  of  these  problems  can  be  obtained  for 
2 

known  a  (since  then  the  "width"  of  the  intervals  is  fixed,  and  not 
a  random  variable).  We  shall  focus  on  one-sided  confidence  intervals 

in  Bechhofer  and  Tamhane  [4],  [5],  and  in  those  papers  shall  give  the 
plans  for  optimal  designs  for  the  MCC  problem. 


T 
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In  the  remainder  of  this  section  we  shall  describe  a  simple  rule 

which  can  be  used  to  determine  whether  or  not  one  design  is  inadmissible 

w.r.t.  another  design  for  given  (p,k,b),  in  the  sense  described  below. 

To  illustrate  the  concept  of  inadmissibility  we  consider  the  case 
.  .  2 

of  one-sided  confidence  intervals;  a  is  assumed  to  be  known.  We 
limit  consideration  to  confidence  intervals  of  the  form 
{aQ-a.  >  a^-a.  i  d  (l<.i<p)}  where  d  >  0  is  a  specified  "yardstick" 
associated  with  the  common  width  of  the  confidence  intervals.  For 
given  (p,k,b)  we  say  that  BTIB  design  is  inadmissible  with 

respect  to  3TI3  design  if  yields  a  larger  joint  confidence 

coefficient  than  does  D0  for  eve ry  d  and  a-  The  following  theorem 
gives  a  characterization  of  inadmissibility  which  is  easy  to  verify. 

Theorem  5.1:  For  given  (p,k,b)  consider  two  BTIB  designs  and 

2  2 

with  parameters  (n^P-^)  and  (n25p0)»  respectively.  Design 

.  .  .  2  2 
is  inadmissible  w.r.t.  design  if  and  only  if  and  >  p^ 

with  at  least  one  inequality  strict. 

Proof  of  sufficiency:  The  probability  associated  with  the  joint  confidence 
statement  {aQ-cu  >_  aQ-cu-d  (l<i<p)}  for  a  design  with  parameters  (n  ,p) 
can  be  written  as 


P{Z^  <  d*N/on  (l<i<p)} 


(5.1) 


where  (Z^,...,Z  )  has  a  p-variate  equicorrelated  standard  normal  distri¬ 
bution  with  common  correlation  p.  As  n  decreases  for  fixed  d  and  p 
(p  increases  for  fixed  d  and  n)  the  probability  (5.1)  increases. 

(The  monotonicity  w.r.t.  p  follows  from  Slepian's  inequality.) 
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Proof  of  necessity:  Suppose  that  the  confidence  coefficient  associated 

with  is  larger  than  that  associated  with  for  every  d  and  a. 

2  2 

Then  from  letting  d  t  ®  (d  +  0).  QED 


If  a  design  is  not  inadmissible  then  it  is  said  to  be  admissible. 

For  any  given  (p,k,b)  one  would  then  limit  consideration  to  all  admissible 
designs.  As  indicated  earlier,  the  problem  of  enumerating  this  set  of 
admissible  designs  for  general  (p,k,b)  appears  to  be  a  formidable  one. 

To  illustrate  the  concept  of  admissibility  we  consider  the  case 

(p,k,b)  =  (3,3,3)  for  which  Dg  and  of  (3.6)  are  two  of  the  pos- 

2 

sible  3TI3  designs;  DQ  has  (n0,pQ)  =  (81/10,  1/3)  while  D1  has 
2 

(n1  ,p1  )  =  (27/2,0).  Thus  is  inadmissible  relative  to  D  .  Actually 

Dq  is  the  only  admissible  design  for  b  =  3,  and  hence  it  is  optimal 
for  b  =  3. 


As  another  example  we  consider  the  case  (p,k,b)  =  (3,2,9)  for 
which  the  only  possible  3TI3  designs  are 


000000000 


Do  = 


(5.2a) 


123123123 


112000000 


D1  = 


(5.2b) 


233123123 


112112000 


D  = 
2 


(5.2c) 


233233123 


In  (5.2),  Dq  has  (ng,pQ)  =  (12,0),  D1  has  (n^p^  =  (108/10,  1/3), 
ar.d  Dj  has  (108/7,  2/3).  Thus  Dg  is  inadmissible  being  dominated 


by  and  are  both  admissible. 


1 


i  m  u  .nttrTTfc  ii  iJfl 
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The  admissibility  criterion  described  in  the  paragraphs  above 
applies  equally  well  to  the  case  of  joint  two-sided  confidence  inter¬ 
vals;  then  the  monotonicity  w.r.t.  p  follows  from  Sidak's  [23] 
results.  (Of  course,  the  optimal  designs  might  be  different  in  the 

one-sided  and  two-sided  cases  for  the  same  (p,k,b)  and  d. )  The 

2 

same  general  ideas  carry  over  for  unknown  o  except  that  then  one 
would  have  to  specify  the  expected  common  "width"  of  the  confidence 
intervals . 

2 

Remark  5.1:  We  note  from  (4.3)  and  (4.4)  that  n  is  a  decreasing 
function  of  both  A^  and  A^  while  p  is  a  decreasing  function  of 
Ag  but  an  increasing  function  of  A  .  Thus  if  for  given  (p,k,b) 
two  designs  have  the  same  values  of  A^  but  different  values  of  A^, 
then  the  design  with  the  smaller  A^  is  inadmissible.  This  phenom¬ 
enon  is  illustrated  by  the  inadmissibility  of  in  (3.6). 

5 . 2  Equivalent  designs 

2  2 

If  for  given  (p,k,b)  two  designs  and  have  ( ,o^)  =  ( n2  ,P? ) , 

then  we  say  that  D  and  D0  are  equivalent .  It  is  easy  to  verify 

that  a  necessary  and  sufficient  condition  for  and  D„  to  be  equiv¬ 

alent  is  that  (A^^A^1^)  =  (Aq^,A^^)  where  (A^\\^^)  refers 
to  design  ( i  =  1 , 2 ) .  A  trivial  case  of  equivalent  designs  arises 

when  one  design  is  obtained  from  another  by  simply  permuting  the  test 
treatment  labels.  In  such  a  case  and  will  be  referred  to  as 

isomorphic  designs;  this  case  is  of  no  practical  interest,  and  we  do  not 
consider  it  further.  A  more  interesting  case  arises  when  the  two  designs 
are  not  isomorphic  and  yet  are  equivalent;  in  our  papers  the  term  equiv¬ 
alent  will  be  reserved  only  for  the  latter  designs. 
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As  an  illustration  of  such  designs  for  (p,k,b)  =  (4,3,7)  consider 
designs  of  (3.3)  and  (say)  of  (3.10d).  The  designs  D9  and 

P,  are  equivalent  since  (A^,X^)  =  (2,2)  for  both  designs;  but  the  de¬ 
designs  are  not  isomorphic.  Equivalent  designs  can  sometimes  provide  flex¬ 
ibility  to  the  experimenter  without  changing  the  confidence  coefficient. 

For  example,  the  experimenter  might  prefer  to  ^  ^  the  con_ 

trol  treatments  are  more  readily  available  for  experimentation  than  any 
of  the  test  treatments. 

As  another  example  of  two  equivalent  designs  which  are  not  isomorphic, 
for  (p,k,b)  =  (6,3,35)  consider 


D 


1 


r 

0  0 

5  replicates  of  /  1  2 

3  6 


0  12  3  1 

4  2  3  4  5  ) 

5  4  5  6  6 


(5.3) 


=  (all  distinct 


<?3> 


combinations  of  (0,1,. ..,6)}, 


(5.4) 


5 . 3  Relationship  to  other  optimality  criteria 

Kiefer  (along  with  many  others)  has  studied  extensively  the  problem 
of  optimal  design  in  a  series  of  papers  starting  with  the  one  [16]  in 
1958.  The  three  main  criteria  considered  by  Kiefer  are  D-,  A-  and  E- 
optimality  which  correspond,  respectively,  to  minimizing  the  determinant, 
trace,  and  the  maximum  eigen  value  of  the  variance-covariance  matrix  of 
the  BLUE  of  the  parameter  vector  of  interest  which  in  our  case  is 
(aQ-a^ , . . . ,a^-a^ ) ' .  The  eigen  values  (neglecting  constant  proportionality 
factors)  of  the  variance-covariance  matrix  of  (a  -ot  , . . . ,a  -a  )' 
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can  be  obtained  from  (4.2)-(4.4);  they  are  (5lg+P*1  )~~  of  multiplicity 
d-1,  and  X  ^  of  multiplicity  1.  Thus  for  given  (p,k,b)  the  three 
criteria  can  be  stated  as  a)  D-optimality :  minimize  {X^CXq+pX^  )P_"}  *, 
i.e.,  maximize  { X Q C A 0-*-pA ^ )*  };  b)  A-optimality : 

minimize{ +  (p-l)( XQ+pXi )  c)  E-optimality:  minimize  X  \  i.e., 

maximize  X  . 

Although  these  criteria  are  simpler  than  ours  it  should  be  kept  in 
mind  that  they  refer  to  an  ellipsoidal  joint  confidence  region  for 
(oIq-o^  , . . .  ) ' .  In  the  case  of  MCC  such  a  confidence  region  is  of 

less  interest  than  the  "rectangular"  confidence  regions  that  we  have  pro¬ 
posed,  and  are  commonly  used.  Hence  we  do  not  consider  the  Kiefer  criter¬ 
ia  further  in  this  series  of  papers. 
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6.  CONCLUDING  REMARKS  AND  DIRECTIONS  OF  FUTURE  RESEARCH 
In  the  present  paper  we  have  introduced  a  new  general  class  of  incom¬ 
plete  block  designs  which  are  appropriate  for  use  in  the  MCC  problem.  We 
refer  to  these  as  balanced  treatment  incomplete  block  (BTIB)  designs.  The 
basic  results  concerning  the  structure  of  such  designs  are  derived,  and 
the  properties  of  the  relevant  estimates  obtained  with  such  designs  are 
given.  Admissibility  and  inadmissibility  of  these  designs  are  defined,  and 
these  criteria  are  used  to  eliminate  inferior  designs. 

The  combinatorial  problem  of  constructing  all  BTIB  designs  for  given 
(p,k,b),  and  the  procedure  for  choosing  an  optimal  design  from  such  a 
set  are  not  solved  in  the  present  paper.  However,  some  methods  of  design 
construction  are  given.  The  aforementioned  problems  are  related  in  the 
sense  that  to  solve  the  optimization  part  completely  one  must  have  con¬ 
structed  most  if  not  all  generator  designs  for  given  (p,k);  the  problem 


of  determining  how  many  generator  designs  exist  for  arbitrary  (p,k), 
and  then  of  enumerating  them  appears  to  be  a  very  formidable  one.  Al¬ 
ternatively,  the  problem  cf  construction  of  BTIB  designs  for  arbitrary 
(p,k,b)  can  be  set  up  in  the  manner  of  Foody  and  Hedayat  [11,  Lemma 
4.1]  which  is  suitable  for  a  solution  on  a  computer.  However,  there 
is  no  guarantee  that  all  BTIB  designs  can  be  generated  in  this  way. 

Also,  the  magnitude  of  difficulty  of  our  problem  is  substantially  great¬ 
er  than  theirs,  and  therefore  a  solution  via  this  route  looks  rather 
remote  at  this  stage. 

For  p  >  2,  k  =  2  and  p  =  3,  k  =  3  the  situation  is  very  simple 
since  it  is  necessary  to  consider  essentially  only  two  generator  designs 
for  each  (p,k)  case;  these  cases  are  considered  in  detail  in  Bechhofer 
and  Tamhane  [4],  and  the  optimal  design  is  given  for  a  large  range  of 
the  useful  (p,k,b)-  and  d-values.  For  p  =  4,5,6  and  k  =  3,4  it 
is  possible  to  enumerate  most  of  the  generator  designs;  these  cases  will 
be  considered  in  Bechhofer  and  Tamhane  [5]  along  with  the  associated  op¬ 
timal  designs. 
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APPENDIX 


Proof  of  Theorem  3.1.  For  given  (p,k,b)  consider  an  arbitrary  design 


with  the  incidence  matrix  {r. J.  Then  it  is  well  known  (see  equation 

(3.1)  of  [16])  that  the  information  matrix  =  (m.  .  }  of 

1112 

$  =  (',o,ai,'"*otp)'  is  given  by 


m. 

i 


1X2 


1  V  2  /•  .  . 

-  —  I  r.  .  (l  =i  ) 

*1  k  j=i  V  1  2 


V2 


(W' 


(A.l) 


Note  that  M:  (p+l)x(p+l)  is  a  singular  matrix  and 

each  i^.  We  require  the  information  matrix  of 

Uct  =  (a  -a  -a  )'  where  U  =  (u.  .  }:  D  x  ( p+1 ) 

■vo  0  1  °P  ^  ‘ 


i2=° 


m.  .  =  0 

1112 


is  given  by 


for 


u. 

l 


112 


1  (i2=1’ 

-1  (i2-ij+^’  •>  •  •  •  >P ) 

0  otherwise. 


In  order  to  avoid  computing  the  generalized  inverse  of  M  we  shall  fol¬ 
low  the  method  given  below  to  obtain  the  desired  information  matrix.  Let 
Q:  p  x (p+1)  be  any  matrix  the  rows  of  which  form  p  orthogonal  con¬ 
trasts.  Then  we  can  write  as  U  =  ^  where  P:  p  x  p  is  a  nonsingular 

matrix.  The  information  matrix  of  Qa  is  given  by  QMQ'.  Therefore  the 

2 

variance-covariance  matrix  (except  for  the  common  factor  a  )  of  Oa 

w 

is  given  by  •  Hence  the  information  matrix  M*  of  Ua  is 


where  V  =  Q'P 


matrix.  It  can  be  easily  verified  that  V 


p  is  given 


Substituting  for  V  in  (A. 2),  the  information  matrix 


can  be  written  as 


where  (A. 3)  follows  from  the  fact  that  the  rows  and  columns  of  M  sum 


to  zero.  For  a  design  to  be  BTIB,  the  matrix 


must  be  completely  syn 


metric,  i.e.,  all  diagonal  elements  of  M*  must  be  equal,  and  all  off 
diagonal  elements  must  be  equal.  Therefore  we  have 


25 


m 


11 


m 


22 


m 

PP 


and 


(A. 4) 


m 


12 


m 


13 


m 


P-l.P- 


Using  expression  (A.l)  for  m.  .  in  (A. 4)  implies  (3.1).  Denoting 

1112 

the  common  value  of  the  XQ^  (l^i^p)  by  X  ,  and  the  common  value 

of  the  X  ^  (i^i^;  by  X1  WS  See  that 

m.  .  =  -X  /k  (i  ^i  ;  l^i-.i  <p)  and  m. .  =  { X  +  (p-1  )X .  }/k  (l<i<p). 

l^l0  1  12  1  2^  11  0  1  = 

Thus 

M*  =  {(X  +pX,  )I  -  X  J}/k 
a.  0  1  l'v 


where  p x p  is  a  matrix  consisting  only  of  l’s.  The  inverse  of 
IJ*  is  the  desired  variance-covariance  matrix  of  Ua;  the  expressions 

for  (4.2)  and  (4.4)  are  then  easily  obtained  from  (M*) 

*\/ 

Derivation  of  the  expressions  for  the  estimates 

The  normal  equations  for  the  least  squares  estimates  (BLUE) 
li,a^,8j  of  u  ,a^ , 8 j  ,  respectively,  (CKi<p,  l^j<b)  are: 

P  b 

N  y  +  7  r.a.  +  k  7  B .  =  G,  (A. 6) 

i=0  1  1  j=l  ^ 


b 

r.y  +  r.a.  +  7  r..B.  =  T.  (0<i<p) 

i  i  i  ij  ]  i  =  - 


P 

ky  +  7  r..a.  +  kg.  =  B.  (l<i<b) 

i=0  ^  1  3  3 


(A. 7) 


(A. 8) 
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where  T.  and  B.  are  defined  in  Section  4.1,  and  G  = 
Substituting 


l  T.  =  l  B.. 
i=0  1  j=l  3 


.  .  P 

U+B.  =(B.  -  T  r..a. )/k 

3  ’  i=0  13  1 


(A. 9) 


from  (A. 8)  into  (A. 7)  we  obtain 


? .  a .  +  71  l  r .  .  B .  -  r-  /  a  f  r,  .  r .  .  =  T .  . 
ii  k  i:  :  k  h^0  h  ^  h:  13  i 


(A. 10) 


Now  I  r..B.  =  B!;  and  for  h  i  i  we  have  that  \  r,.r.  .  =  A  (A  ) 

j=l  ^  3  1  j=l  h)  01 

if  only  one  of  h  and  i  =  0  (if  both  h  and  i  i  0).  From  (A. 10)  for 


i  =  0  we  obtain 


g  J*  ^  ^ 

Vo  +  T  "  £{“o  J  roj  +  xo  V  =  V 

j=l  J  h=l 


(A. 11) 


P  .  i  b  2 

Using  the  fact  that  T  a,  =  -  a.  and  that  r  -  —  j  r . 
&  ,  ,  h  0  0  k  ,L.  0i 

h=l  3=1  J 

from  (A. 11)  we  obtain 


m _ =  pA./k, 

00  ^  0 


a,  = 


0  (p+l)A  ' 


(A. 12a) 


In  the  same  way  we  obtain 


V'Wo  . 

“i  -  (liia,)' 


(A. 12b) 


Combining  (A. 12a)  and  (A. 12b)  we  obtain  (4.1). 
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Derivation  of  the  formula  for  the  adjusted  treatment  sum  of  squares 
in  the  analysis  of  variance  table  (Table  4.1) 

Following  the  Scheffe  [22]  notation,  let  denote  the  minimum 

error  sum  of  squares  (SS)  under  the  assumptions  £2  of  Section  2,  and 
let  S  denote  the  minimum  error  SS  under  ui  =  H  n  i]  where  the 

U) 

hypothesis  is  H:  =  0  (O^iip).  Then  the  SS  of  treatments  adjusted 

for  blocks  is  given  by 


SS  w  .  =  S  -  S^. 
treat ( ad j.)  u  ft 


(A. 13) 


Sn  =  SS 
ft  error 


d  b  k 

=  I  ll  (y 


i=0  j=l  h=l 


ijh  i  j  ijh 


(A. 14) 


where  =  1  if  the  ith  treatment  is  assigned  to  the  hth_  plot  of  the 

jth_  block  (04i<_p,  l£j<£,  Tlh^k)  and  -  0  otherwise.  Substituting  6. 
from  (A. 9)  into  (A. 14),  and  expanding  and  simplifying  the  resulting  square 
we  obtain 


S  =  H  -  i  l  B2.  +  \  r.a? 

ft  k  ,L.  i  .Ln  i  i 

j=l  J  i=0 


■r  l  (I  r..0L.V  .  Q 

k  j=l  i=0  1 


(A. 15) 


P  b  k  2  2  ?  . 

where  H  =  Y  Y  Y  and  Q  =  —  T  o.Q.. 

i=oj=i  hu  i5t; ijh  k  iio  11 

U  =  G/N  and  (L  =  B^/k  -  u;  hence 


Under  ui  we  have 


r 
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Soj  =  H  -  i  l  B2.. 

k  j=l  j 


(A. 16) 


Subtracting  (A. 15)  from  (A. 16)  we  have  from  (A. 13)  that 


SS 


treat ( ad i ) 


•>i  =  Q  +  r  I  (  1  r..a.)2  -  \  r.a? 
''  k  j‘l  i=0  «  1  i=0  11 


2  lr2"2  2?  ?  * 

=  Q  -  I  (r.  --  £  r..)a.  +  —  £  £  a.  a.  A  . 

i=0  1  k  ^  ^  1  k  <  =n  <  =•  *i  ^  lo  M 


j=l 


il=°  i2=ii+1  1  212 


pA 


=  Q 


0  *2  <V(p-1)Al}  ?  -2 
~  ao - Z -  l  «■ 


k  0 


i=l 


(A. 17) 


2A, 


2A  a  2 

0  0  v  1 

+  — - )  a.  +  — —  /  )  a.  a. 

k  .  ,  l  k  .  ,  .  r  ..  l,  i„ 

i=l  l  =1  i2=i1+l  1  2 


<Pt2>V0  VV  ?  -2  X1  ,  ?  -  ,2 

Q - k - — k -  J,  “i  +  T  (J,  V 

1=1  1=1 


L(p+2)AQ-A1]a2  (AQ+pA1) 

=  Q  _  _ 

1  ^  O 

In  the  above  we  have  used  the  relations  r.  -  —  /  r? .  =  m. .  (from  (A .1)) 

i  k  i]  ii 

where  mQ0  =  pA0/k,  m.^  =  { \Q  +  (p-1  )A^ }/k  (l<><p),  and  ^  a.  =  -otQ. 

i=l 

Substituting  in  (A. 17)  for  from  (A. 12a)  and  for  a.  (l<iip)  from 

(A. 12b),  and  after  some  tedious  algebra  we  obtain  the  expression  for 

SS  ....  as  given  in  Table  4.1.  The  other  SS  expressions  and  the 

treat(ad') 

df  in  the  table  are  obtained  in  a  straightforward  way. 
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